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We call monomer a B-DNA base-pair and examine, analytically and numerically, electron or hole 
oscillations in monomer- and dimer-polymers, i.e., periodic sequences with repetition unit made of 
one or two monomers. We employ a tight-binding (TB) approach at the base-pair level to readily 
determine the spatiotemporal evolution of a single extra carrier along a N base-pair polymer. We 
study HOMO and LUMO eigenspectra as well as the mean over time probabilities to find the carrier 
at a particular monomer. We use the pure mean transfer rate k to evaluate the easiness of charge 
transfer. The inverse decay length fi for exponential fits k{d), where d is the charge transfer distance, 
and the exponent rj for power law fits k{N) are computed; generally power law fits are better. We 
illustrate that increasing the number of different parameters involved in the TB description, the 
fall of fc(d) or k{N) becomes steeper and show the range covered by /3 and rj. Finally, both for the 
time-independent and the time-dependent problem, we analyze the palindromicity and the degree of 
eigenspectrum dependence of the probabilities to hnd the carrier at a particular monomer. 


I. INTRODUCTION 

Recently, we studied B-DNA dimers, trimers and poly¬ 
mers (we call monomer a B-DNA base-pair), with a tight- 
binding (TB) approach at the base-pair level, using the 
relevant on-site energies of the base-pairs and the hop¬ 
ping parameters between successive base-pairsi*^. Our 
method allows us to readily determine the spatiotempo¬ 
ral evolution of holes or electrons along a N base-pair 
DNA segment through the solution of N coupled dif¬ 
ferential equationsii^. We showed that for all dimers 
and for trimers made of identical monomers the carrier 
movement is periodic with frequencies in the mid- and 
far-infrared i.e. approximately in the THz domainii^, a 
region of intense research^. Increasing the number of 
monomers above three, periodicity is generally lostii^. 
Even for the simplest tetramer, the carrier movement is 
not periodic^. For periodic cases, we defined^i^ the max¬ 
imum transfer percentage p (e.g. the maximum prob¬ 
ability to find the carrier at the last monomer having 
placed it initially at the first monomer) and the pure 
maximum transfer rate ^ = pf {T being the period, 
/ the frequency). For all cases, either periodic or not, 
the pure mean transfer rate k (cf. Eq. ITT]) and the speed 
u = kd, where d = {N — 1) x 3.4 A is the charge transfer 
distance, can be used to characterize the system. Our 
analytical calculations and numerical results show that 
for dimers k = 2^ and for trimers made of identical 
monomers k « 1.3108p. Using k to evaluate the easiness 
of charge transfer, one can calculate the inverse decay 
length P for exponential fits k{d) and the exponent p for 
power law fits k{N). Studying a few polymers and seg¬ 
ments taken from experiments^, we found that /3 ~ 0.2 - 
2 and p « 1.7 - 17. 

Carrier oscillations within “molecular” systems have 
been occasionally studied in the literature. Real- 
Time Time-Dependent Density Functional Theory (RT- 
TDDFT)^ simulations predicted oscillations (~ 0.1-10 
PHz) within p-nitroaniline and FTC chromophore^, zinc 


porphyrin, green fluorescent protein chromophores and 
adenine-thymine base-pair—. In a simplified single- 
stranded helix of 101 bases, subjecting the system to 
a collinear uniform electric field, THz Bloch oscillations 
could be induced^. Single and multiple charge transfer 
within a typical DNA dimer in connection to a bosonic 
bath, where each base-pair is approximated by a single 
site, as in our TB approach, has been studied in Ref.—. In 
the subspace of single charge transfer between base-pairs 
and having initially placed the charge at the donor site, 
the authors obtain a period slightly greater than 10 fs, 
having used a “typical hopping matrix element” 0.2 eV. 

Using our equationl^ / = A = putting t = 

0.2 eV for the “typical hopping matrix element” and iden¬ 
tical dimers i.e. difference of the on-site energies A = 0, 
we obtain a period T Ri 10.34 fs in accordance with the 
dotted line in Fig. 4 of Ref—. 

In the present article, we focus on periodic DNA poly¬ 
mers with a repetition unit made of one or two monomers 
and analyze the parameters which favor charge transfer. 
A synopsis of the theory behind the calculations is given 
in Section El both for the time-dependent and the time- 
independent problem. In Section Hill we distinguish three 
special types of polymers, we discuss analytical solutions, 
and we present and discuss our numerical and analytical 
results. Both for the time-independent and the time- 
dependent problem, we introduce two important proper¬ 
ties of the probabilities to find the carrier at a particu¬ 
lar monomer: palindromicity and degree of eigenspectrum 
dependence. Moreover, we use k to characterize the eas¬ 
iness of charge transfer and compare exponential fits of 
k{d) with power law fits of k{N). We also compare hts in¬ 
cluding any number of monomers with fits including only 
odd or only even numbers. In Section llVI we summarize 
our conclusions. Finally, we just mention that although 
for A > 4 the oscillations are not periodic, Fast Fourier 
Transform analysis shows that the frequency content re¬ 
mains strong in the THz domain. 
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II. THEORY 

By YX we denote two successive base-pairs, according 
to the convention 

5' 3' 


^'h/l' HOMO/LUMO on-site energy of base-pair 

fjL, and ^ is the hopping parameter between base- 

pair fi and base-pair fi'. The values of and ^ 

used in the present work are the same employed in 

RefsJ-i^. 

To solve Eq. |T] we define the vector matrix 


Y 

X 

3' 


Ycompl 

Xcoinpl 


5' 


: ( 1 ) 

for the DNA strands orientation. We denote by X, 
Xcoinpl; X, Ycompi DNA bases, where Xcompi (Xcompi) is 
the complementary base of X (Y). In other words, the 
notation YX means that the bases Y and X of two suc¬ 
cessive base-pairs are located at the same strand in the 
direction 5' — 3'. X-Xcompi is the one base-pair and Y- 
Ycompi is the other base-pair, separated and twisted by 
3.4 A and 36°, respectively, relatively to the first base- 
pair. For example, the notation GT denotes that one 
strand contains G and T in the direction 5' — 3' and the 
complementary strand contains C and A in the direction 
3' - 5'. 


For each base-pair or monomer, the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) play a key role, since we sup¬ 
pose that an extra hole or electron inserted in a DNA 
segment travels through HOMOs or LUMOs. 


x{t) 


Then, Eq. breads 


Ai{t) 

Mt) 

A^it) 


x{t) = Ax{t), 


(5) 


( 6 ) 


where A = —-^A and the matrix A is a symmetric tridi¬ 
agonal matrix shown in the Appendix We solve Eq. [6] 
using the eigenvalue method, i.e. looking for solutions of 
the form x{t) = ve^* x{t) = Xve^*. Hence, Eq.jSlreads 


Av = Xv 

or, with A = —f A, 

Av = Xv, 


(7) 

( 8 ) 


i.e. we have to solve an eigenvalue problem. Provided 
that the normalized eigenvectors iTk corresponding to the 
eigenvalues A^ of Eq. [5] are linearly independent (which 
holds in all cases studied), the solution to our problem is 


A. Time-dependent problem 

For a TB description at the base-pair level, the time- 
dependent single carrier (hole/electron) wave function of 
the DNA polymer, 4'2y^(r, t), is written as a linear com¬ 
bination of base-pair wave functions with time-dependent 
coefficients, i.e.. 




N 

x{t) = 

k=l 


(9) 


If we initially place the carrier at base-pair 1 and we 
want to see how the carrier will evolve, time passing, then 
the initial condition would be 


( 10 ) 





■ 7li(0) ■ 


■ 1 ■ 

N 

Y,AAt) <«(r). 

(2) 

f(0) = 

A2(0) 

= 

0 

fl^l 



_ An{0) _ 


0 


'L^y^^(r) is the base-pair’s HOMO or LUMO wave 
function (H/L). The sum is over all base-pairs of the 
DNA polymer. |A^(t)p gives the probability to find the 
carrier at the base-pair /i, at the time t. 

Using the time-dependent Schrodinger equation 

L \ HI {jDNA,^,DNAt„ fo\ 

- =H ^h/lAA) (3) 

and Eq. [2] as well as the details described by Hawke et 
al.— , we find that the time evolution of the coefficients 
A^(t) obeys the TB system of differential equations 


ih 


dAfj^ 

dt 


= E 


bp(v) 

H/L 


A I a _i_ Apiu-’V+A A 


H/L 


V-l-1' 


(4) 


However, one could initially place the carrier at another 
base-pair (hence 1 would be at the corresponding place at 
the right-hand side) or even imagine to initially distribute 
the carrier probability equally among monomers (hence 
all right-hand side components would be 1/a/N). From 
the initial conditions we determine c^. 

An estimation of the transfer rate can be obtained^ as 
follows: Supposing that initially i.e. for t = 0 we place 
the carrier at the first monomer lEg. fTUl) . then |Ai(0)p = 
1, while all other \Aj (0)p = 0, j = 2,..., N. Hence, for a 
polymer consisting of N monomers, a pure mean transfer 
rate can be defined as 


{\ANm 

tN mean 


( 11 ) 
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where tMmean is the first time |^ 7 v(t)P becomes equal 
to (|A 7 v(t)P) i.e. “the mean transfer time”. Finally, the 
speed of charge transfer could be defined as u = kd; d = 
(TV — l)x 3.4 A is the charge transfer distance. 


B. Time-independent problem 


The time-independent Schrodinger equation 

HDNA^DNA^r) = Avl/gf/(r) (12) 

can be solved expanding the time-independent single car¬ 
rier (hole/electron) wave function of the DNA polymer, 
as a linear combination of base-pair wave func¬ 
tions with time-independent coefficients, i.e.. 


N 

= (13) 

fj.=l 

I P gives the probability to find the carrier at the base- 
pair fj,. The problem specified in Eqs. (ini-dTSl) is equiv¬ 
alent with Av = \v (Eq. [5|) , with 


hi 

Fa 


Tat 


(14) 


In other words, the eigenvalues and eigenvectors of Eq.fT^ 
are Afc and Vk, respectively; v^k = 


so, we only have one non-zero TB parameter. In type /?', 
still A = 0 but A so, we have two non-zero TB 
parameters. Finally, in type 7 ', A A 0 and A so, 
we have three non-zero TB parameters. 

The eigenproblems we have to solve refer to a tridiag¬ 
onal Toeplitz matrix of order N for type a' polymers (cf. 
Eq. 1251 the analytical solution is rather simple) and a 
tridiagonal 2-Toeplitz matrix of order N for type /3' (cf. 
Eq. [31]) and 7 ' (cf. Eq. |35|) polymers. These eigenprob¬ 
lems have been studied in RefAi where the characteristic 
polynomial of a tridiagonal 2-Toeplitz matrix is shown to 
be closely connected to polynomials satisfying the three 
point Chebyshev recurrence formula - an extension of the 
well-known result for a tridiagonal Toeplitz matrix. Two 
theorems (2.3 and 2.4) describe the eigenvalues for odd 
and even N—. When N is odd the eigenvalues can be 
expressed explicitly in terms of Chebyshev zeros^i. Al¬ 
though for even N there is no explicit formula, a recipe 
to produce the eigenvalues is givenii. Specifically, these 
theorems refer to the tridiagonal 2-Toeplitz matrix of or¬ 
der n, given by Eq. 2.8 of RefAi (for us n = N.): 

ai Pi 

7i 0^2 P 2 

0 72 ai 

— 

71 

0 

Theorem 2.3 of Ref The eigenvalues of the tridiagonal 
2-Toeplitz matrix of order 2 to-|- 1 given in Eq.|TS](Eq. 2.8 
of Ref.—) are ai and the solutions of the quadratic equa¬ 
tions 


Pi 

Q !2 


(15) 


III. PERIODIC POLYMERS WITH 
REPETITION UNIT MADE OF ONE OR TWO 
MONOMERS 

Let us focus on periodic DNA polymers with repetition 
unit made of either one monomer or two monomers (a 
dimer). We distinguish three types of DNA polymers: 
(type a') poly(dG)-poly(dC) and poly(dA)-poly(dT), 
(type P') GCGC..., GGCG..., ATAT..., TATA..., and 
(type 7') TCTC... = GAGA..., CTCT... = AGAG..., 
AGAG... = GTGT..., GAGA... = TGTG... . 

Let us define A := - E'’p^‘='> \, where E^p<^°'> is the 

on-site energy of the carrier at odd monomers (^=1,3, 
5, ...) and is the on-site energy of the carrier at 

even monomers (/i = 2, 4, 6 , ...). Let us by the way define 
E := E^P^°'^ + E^P'^^P Further, counting from the start, 
let us call t^P the hopping parameter from odd to even 
monomers (between ^ = 1 and /r = 2 ...) and P’p the 
hopping parameter from even to odd monomers (between 
fjL = 2 and /x = 3 ...) . For simplicity, we have dropped 
the indices H/L. 

Then, we realize that the intricacy of the energy struc¬ 
ture - i.e. the number of different parameters involved 
in the TB description - increases from type a' to type P' 
and further to type 7': In type a', A = 0 and t^P = P’p, 


(oi - A) (02 - A) - / 3 i 7 i -I- V/ 3 i/ 327 i 72 T’r + P2I2 = 0, 

(16) 

where = 2 cos r = 1 , 2 ,..., m, are the zeros of 
Pmid) defined by Equations IT^ and EUl 
Theorem 2.4 of Ref The eigenvalues of the tridiagonal 
2-Toeplitz matrix of order 2m given in Eq. jTSKEq. 2.8 of 
Ref—) are the solutions of the quadratic equations 


(ai - A)(02 - A) - Pi^i + PlP2lll2Qr + P2I2 = 0 , 

(17) 

where Qr, r = 1 , 2 ,..., m, are the zeros of q'm[d) defined 
by Equations 1121 and HU 

Equations [TSlfTOl are the Chebyshev three point recur¬ 
rence formula. 

Pm+l{d) = dPm{T)-p'm-l{d), ( 18 ) 

Wi(m) = m'mid) - d'm-lid)- (19) 

The initial polynomials are 

Po(Ai) = l ( 20 ) 

9o(m) = 1 q'iip) = T + P- ( 21 ) 
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Finally, 


^2 ^ ^272 
^ Pill 


( 22 ) 


M = 


p 


(23) 


(ai — A)(a 2 — A) 
/3i7i 


(24) 


The eigenvectors are found in terms of polynomials sat¬ 
isfying the three point recurrence relationship^^. For 
the tridiagonal 2-Toeplitz matrices of our interest, up 
to our knowledge, explicit eigenvalues have been found 
only for odd TV— which agree with the results of Ref—. 
Throughout this work we solve the eigenproblems numer¬ 
ically and additionally we compare to analytical results. 


A. Stationary states (time-independent problem): 

Eigenvalues and Eigenvectors 

In this work we calculate the eigenvalues and eigenvec¬ 
tors numerically. However, in some cases, we compare 
with analytical solutions. To simplify the notation, be¬ 
low we write E^p = E, P’p = t and P'p = t'. 


1. type a' (and type a' cyclic) 

For type a' [poly(dG)-poly(dC) and poly(dA)- 
poly(dT)], the matrix A is a symmetric tridiagonal uni¬ 
form matrix 


r ^ 

t 

0 •• 

• 0 

0 

0 1 

t 

E 

t • • 

• 0 

0 

0 

0 

0 

0 •• 

■ t 

E 

t 

L 0 

0 

0 •• 

■ 0 

t 

e\ 


with eigenvalues 

Afc =E + 2tcos . (26) 

where k = 1,2,...,iV. All eigenvalues are real and 
distinct (non degenerate) since the matrix is symmet¬ 
ric (A = a"*"), all eigenvalues are symmetric around A, 
for odd N the trivial eigenvalue (= E) exists, and all 
eigenvalues lie in the interval {E — 2t,E + 2t). The eigen- 
spectrum of type a' polymers is shown in Fig. [T] The /i 
component of the k eigenvector is given by 


— 


r^~ . 
V TV + 1 


pLkTT A 

TV + lj ’ 


(27) 


where k = 1,2,...,TV and /r = 1,2,...,TV. Since Vf^k 
do not depend on E or t, then, for any k, the probabil¬ 
ity to find the carrier at a particular monomer /r, 
also does not depend on E or t. This property (let’s call 
it eigenspectrum independence of the probabilities) 
is conserved in the time-dependent case (cf. subsubsec- 
tion IIII BT|) . Since sin( ) = ±sin(-^^) it fol¬ 

lows that for each eigenstate fc, Ir'/ifcP are palindromes 
i.e. the occupation probability for the p-th monomer is 
equal to the occupation probability of the (TV — p + l)-th 
monomer. This property (let’s call it palindromicity) 
is conserved in the time-dependent case (cf. subsubsec¬ 
tion nr B II) . The eigenspectra of type a' polymers are 
shown in the first two rows of Fig. [T] 

Furthermore, one could imagine the cyclic polymers 
with A(1,TV) = = A(TV, 1) = ^ 0. For 

cyclic type a' polymers, the matrix A is a symmetric 
tridiagonal uniform matrix with two “perturbed corners” 


r ^ 

t 

0 •• 

• 0 

0 

t 

t 

E 

t ■■ 

• 0 

0 

0 

0 

0 

0 •• 

■ t 

E 

t 

. t 

0 

0 •• 

■ 0 

t 

e\ 


whose eigenvalues are 

( 2 \ 


where fc = l,2,...,TV. Generally, all eigenvalues are not 
distinct but degeneracies exist. The number of discrete 
eigenvalues M = for TV odd and M = for TV 
even. Eigenvalues of this and other tridiagonal Toeplitz 
matrices with four perturbed corners can also be found 
in Ref—. The p component of the k eigenvector is given 
by 


where p = 1,2, ...,TV. Since for any 

eigenstate fc, the occupation probability is equal for all 
monomers. The eigenspectra of type a' cyclic polymers 
are shown in the last two rows of Fig. [TJ 


2. type P' 

For type P' polymers, the matrix A is 

' E t 0 0 ••• ■ 

t E t' 0 ■■■ 

A= Q t' E t ■■■ 


(31) 


For odd TV, A has the same number of t and t'. 
Hence, for odd TV, its eigenvalues and eigenvectors 
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have some noteworthy properties: For odd N, for the 
same set of parameters the set of eigenval¬ 

ues {Afc} remains the same if we interchange the se¬ 
quence of base-pairs, i.e. {Afe}(XY...)={Afc}(YX...); e.g. 
{Afc} is the same for HOMO GCGCGCG and HOMO 
CGCGGGG. Moreover, for odd N, for the same set of 
parameters the eigenvectors have the proper¬ 

ties |u^fe(XY...)i = |i;(Ar_^+i)fe(YX...)| and |u^fc(XY...)| = 
k/i(Ar-fe-i-i)(XY...)|. For odd N, the eigenvalues can be 
writteni^ as 


Afc 


E + ^/t'^ + t'^ + 2tt' cos 9k , fc = l,...,m 
E — \/+ 2U' cos 0k , fc = TO -I- 1,..., 2m 
E, k = N 

(32) 


Ok 


2k'7z 

N-\-l 

2{k—m)'7T 

iv+i 


k = 1,..., TO 

fc = TO -|- 1,..., 2 to ’ 


(33) 


which are equivalent with the resulting eigenvalues in 

RefJi 


1 

\E, 

and 

{Afc} = 

1 f;± 

y^"+^'"+2tt'C0s(;^) 


where to = and r = 1, 2,... to. Analytical expres¬ 
sions for the eigenvectors, for odd TV, can be found in 
Refji^. For odd TV, it is worth noting that the eigen¬ 
vectors v^k depend on t and t', hence, for any fc, the 
probability to find the carrier at a particular monomer 
/i, |u^fcp also depends on t and t'. Hence, in contrast 
to type a' polymers, now we have partial eigenspec- 
trum dependence of the probabilities i.e. dependence 
on the hopping parameters but not on the on-site en¬ 
ergy. For odd TV, Ir’/^fcp are palindromes only for even /i; 
this property is conserved in the time-dependent case (cf. 
subsubsection Uli B 21) . For even TV, the situation is more 
complicatedii. We have not encountered an analytical 
solution, in the literature, yet. For even TV, A does not 
have the same number of t and t'. For even TV, |u;jfcp are 
palindromes for all /x; this property is conserved in the 
time-dependent case (cf. subsubsection HH B 21) . Hence, 
we have palindromicity for TV even, but for TV odd 
only partial palindromicity. The eigenspectra of type 
P' polymers for odd and even TV are shown in Fig. [21 


3. type 7 ' 


For type 7 ' polymers, the matrix A is 


A = 


E° t 0 0 ••• 

t E^ t' 0 ••• 

0 t' E° t ■■■ 


(35) 


For odd TV, A has the same number of t and t'. For odd 
TV, the eigenvalues can be written^ as 


1 

\ E°, 

and 

(Afcl = 

[ 1 ± ' 

\J (#)" + + 2 «' cos 


(36) 

where to = ^^ ^ and r = 1, 2,... to. This is in accordance 
with Ref.—. For odd TV, analytical expressions for the 
eigenvectors can be found in Ref—. It is worth noting 
that the eigenvectors depend on E°, E^, t and t', 
hence, for any fc, the probability to find the carrier at a 
particular monomer /i, \v^k\'^ also depends on FT®, FT®, t 
and t'. Hence, in contrast to type a' polymers now we 
have eigenspectrum dependence of the probabilities. 
For even TV, the situation is more complicatedii. We have 
not encountered an analytical solution, in the literature, 
yet. For even TV, A does not have the same number of t 
and t'. The eigenspectra of type 7 ' polymers for odd and 
even TV are shown in Fig. [3] 


B. Mean —over time— Probabilities 
(time-dependent problem) 

The behavior of the mean -over time- probabilities to 
find the carrier at base-pair /x, (|A^(t)p), is different in 
types o', 7 '. 


1. type a' (and type a' cyclic) 

For type a' polymers, {\A^(t)p) are palindromes (no¬ 
tice that Vfik for type a' polymers are palindromes, too) 
and they do not depend on the on-site energies and the 
hopping parameters, but only on TV. In other words, 
eigenspectrum independence and palindromicity 
of the probabilities are reflected here from the station¬ 
ary case (cf. subsubsection IHI A 1 1) . 

If we initially place the carrier at the 1st monomer, 
then the mean -over time- probabilities are 

{\A,{t)f) = (lAv(t)l^) = ^(]^>VTV > 2, (37) 

{\A2{t)f) = ... = (|A^_i(t)|2) = VTV > 3. (38) 

In Fig. m we illustrate, (|A^(t)p) for HOMO and 
LUMO poly(dA)-poly(dT) [the figures are identical for 
poly(dG)-poly(dG)] if we initially place the carrier at the 
1st monomer: at the left column for TV = 5 and TV = 17 
and at the right column for TV = 6 and TV = 18. These 
follow Eqs. (I571) - (I35)) i.e. depend only on TV. In Fig. (TUI 
in the Appendix [B] we display other properties of a char¬ 
acteristic polymer of type a' (poly(dA)-poly(dT)) either 
for hole (left column) or electron (right column) transfer. 
Again, we initially place the carrier at the 1st monomer. 
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Generally, for type a' polymers, for initial placement 
of the carrier at a particular monomer, we obtain 2 (n+i) 
additional mean -over time- probability at the monomer 
where the initial placement is made and at the symmetric 
relative to the polymer center monomer. Hence, for N 
odd, for initial placement at the central monomer, that 
central monomer obtains 2 (jv+i) ^idditional mean -over 
time- probability. In other words, if we call if; and x the 
mean -over time- probabilities at the favored and at the 
rest monomers, respectively, then '0 = X + 2 {n+i) 

■0 = X + 2 {n+i) ^ initial placement at the 

central monomer). Since the sum of all the mean -over 
time- probabilities is 1 , we obtain 




3 1 

2(1V+ 1)’ ^ “ iV + 1’ 


(39) 


except for N odd and initial placement at the central 
monomer in which case we obtain ip = X = 

On the contrary, if we imagine to initially distribute 
the carrier probability equally among monomers, then 
we obtain the mean -over time- probabilities (from edge 
to center monomers) as , jv(jv-i-i) ; ■ ■ ■ i while for 

N odd the mean -over time- probability at the central 
monomer is 

Let us now allow the first monomer to interact with 
the last monomer with i.e. for cyclic type a' poly¬ 
mers. Then, for initial placing of the carrier at a particu¬ 
lar monomer, we obtain additional mean -over time- 
probability at the monomer where the initial placement 
is made and at the diametric monomer if it exists (i.e. 
for even N). In other words, if we call ip and x the mean 
-over time- probabilities at the favored and at the rest 
monomers, respectively, then ip = x + ^ ■ Since the sum 
of all the mean -over time- probabilities is 1 , we obtain 

for even iV, (40) 

2iV - 1 N -1 

^ for odd A^. (41) 

This is depicted in Fig. [S] On the contrary, if we imag¬ 
ine to initially distribute the carrier probability equally 
among monomers, this initial equidistribution is con¬ 
served and no mean carrier movement is observed. 


2. type P' 

Let us put the carrier initially at the first monomer. 
For type /?' polymers, (|H^(t)p) do not depend only on 
N in contrast to type a' polymers, i.e., for type /?' poly¬ 
mers eigenspectrum independence of the probabil¬ 
ities does not hold. However, interestingly, for N even, 
(|H^(t)P) are palindromes, while for N odd, this only 
holds for even /i. In other words, we have palindromic- 
ity for N even, but only partial palindromicity for 


N odd. These symmetry properties, can be summarized 
as 


N 

N 


(42) 

even, * = 0 , 1 ,... iV — 1 , or 
odd, i = 1 ,3,... — 2 . 


In Fig. ini we depict (|H^(t)p) for HOMO and LUMO 
GCGC...; at the left column for A^ = 5 and A^ = 17 and 
at the right column for A^ = 6 and N = 18. For GGGG..., 
the hoping parametersii^ are quite different in magnitude 
for holes but rather similar for electrons, i.e. while for 

holes 1^1 = 1^1 = 0 . 2 , for electrons |^| = |^| = 

1.25. This leads to almost disrupted hole transfer when 
the number of repetition units is not integer i.e. for odd 
N. 


3. type 7 ' 

Let us put the carrier initially at the first monomer. 
For type 7 ' polymers, (|H^(t)p) do not depend only on N 
in contrast to type a' polymers, i.e., for type 7 ' polymers 
eigenspectrum independence of the probabilities 
does not hold. In Fig.[7]we depict (|H^(t)p) for HOMO 
and LUMO AC AC...; at the left column for TV = 5 and 
TV = 17 and at the right column for TV = 6 and TV = 18. 
For HOMO ACAC... where incidentally \t^P\ = |, for 

TV odd, (|A^(T)p) are palindromes. 


C. Pure mean transfer rate fits 

Next, we initially place the carrier at the 1st monomer 
and examine the pure mean transfer rates. Specifically: 
In Fig. [5] (correlation coefficients) and Fig. [5] (/3 and 77 ), 
we compare the exponential fit /c = ( 1 st row), 

the exponential k = A + kQe~^'^ (2nd row, A typically 
results tiny), and the power law ht fc = kQN~P (3rd row) 
for type a', P' and 7 ' polymers. These Hts are carried out 
up to TV = 60, i.e., d = 200.6 A, since it is knowi>i^“— 
that a carrier can migrate along DNA over 200 A. The 
1st column refers to fits including all TV while in the 2nd 
column we fit separately even and odd TV. We observe 
that, generally, the hts improve when we separate even 
and odd TV. 

Furthermore, it is evident that the power law hts are 
signihcantly better. This agrees with the assertio n^^d^ 
that when every single hopping step occurs over the same 
distance, the hopping mechanism is described better by a 
power-law ht. Since here we study periodic polymers, it 
seems that the above scenario holds. For type a', since 
(|A^(t)P), both for even and odd TV, follow the same 
Eqs. (l^ - (l38l) . a htting k{N) does not really depend on 
which TV - even or odd - we include. This is obvious in 
terms of correlation coefficients or rj at the last row of 
Figs. [51 ini This does not hold for types /3' and 7 ' where 
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-since the repetition unit is a dimer- we need to separate 
fits for odd and even N. 

Finally, it is evident that, as a general trend, the fall 
of fc as a function of c? or TV becomes steeper when the 
intricacy of the energy structure is increased i.e. from 
type a' to type /3' and further to type 7 '. 

IV. CONCLUSION 

We have systematically studied electron or hole oscilla¬ 
tions in B-DNA monomer-polymers and dimer-polymers, 
i.e., periodic sequences with repetition unit made of one 
or two monomers, where monomer is a base-pair. We 
used a tight-binding approach at the base-pair level to 
determine the temporal and spatial evolution of a single 
extra carrier along the N base-pair DNA polymer. We 
studied the HOMO and LUMO eigenspectra as well as 
the mean over time probabilities to find the carrier at a 
particular monomer. 

Furthermore, we used the pure mean transfer rate k 
to evaluate the easiness of charge transfer and estimated 
the inverse decay length /3 for exponential fits k{d), where 
d = {N — l)x 3.4 A is the charge transfer distance, and 
the exponent rj for power law fits k{N). It seems that 
the power law fits are significantly better. We also illus¬ 
trated that increasing the energy structure intricacy i.e. 
the number of different parameters involved in the TB 
description, the fall of as a function of d or TV becomes 


steeper, and we showed the range covered by j3 and ry. 

Finally, we combined analytical and numerical solu¬ 
tions for the time-independent and the time-dependent 
problem and analyzed palindromicity and degree of 
eigenspectrum independence of the probabilities to find 
the carrier at a particular monomer. Eigenspectrum inde¬ 
pendence means that the probability to find the carrier at 
a particular monomer does not depend on the on-site en¬ 
ergies and the hopping integrals. Palindromicity means 
that the occupation probability for the /r-th monomer is 
equal to the occupation probability of the (TV — l)-th 
monomer. Type a' polymers display both palindromic¬ 
ity and eigenspectrum independence of the probabilities. 
Type /?' polymers display partial eigenspectrum depen¬ 
dence, and palindromicity for TV even but only partial 
palindromicity for TV odd. Generally, type 7 ' polymers 
do not have either eigenspectrum independence or palin¬ 
dromicity of the probabilities. 
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FIG. 1: Eigenspectra of type a' polymers and type a' cyclic polymers. 
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FIG. 2: Eigenspectra of type P' polymers, 
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HOMO TCTC... 
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FIG. 3: Eigenspectra of type 7 ^ polymers. 
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HOMO poly(dA)-poly(dT) /V=5and/V=17 



LUMO poly(dA)-poly(dT) /V=5and/V=17 



HOMO poly(dA)-poly(dT) A/=6andA/=18 



LUMO poly(dA)-poly(dT) /V=6and/V=18 



FIG. 4: Mean -over time- probabilities (|A^(t)|^) of type oc polymers, if we initially place the carrier at the 1st monomer, -i/; 
and X denote (|7l^(t)p) at one of the favored monomers and at one of the rest monomers, respectively. Typically, we illustrate 
the situation for HOMO and LUMO poly(dA)-poly(dT) but the figures for poly(dG)-poly(dC) are identical. [Left column] 
N = 5 and N — 17. [Right column] N = 6 and N = 18. {[4l^(t)[^) follow Eqs. (I37I) - (I38II i.e. are palindromes and depend only 
on N and not on the on-site energies and the hopping integrals. 


HOMO poly(dG)-poly(dC) cyclic N= 5 and A/= 17 
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HOMO poly(dG)-poly(dC) cyclic N=6 and A/= 18 
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LUMO poly(dG)-poly(dC) cyclic N= 5 and N= 17 
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LUMO poly(dG)-poly(dC) cyclic N= 6 and N= 18 
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FIG. 5: Mean -over time- probabilities ([7l^(t)[^) of type a' cyclic polymers, if we initially place the carrier at the “1st 
monomer”, rjj and x denote {\A^{t)f) at one of the favored monomers and at one of the rest monomers, respectively. Typically, 
we illustrate the situation for HOMO and LUMO poly(dG)-poly(dC) but the figures for poly(dA)-poly(dT) are identical. [Left 
column] N = 5 and N = 17. [Right column] N = 6 and N — 18. {[A^(t)[^) follow Eos. HOETI i.e. depend only on N and not 
on the on-site energies and the hopping integrals. If -on the contrary- we imagine to initially distribute the carrier probability 
equally among monomers, this initial equidistribution is conserved and no mean carrier movement occurs. 
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HOMOGCGC... yV=5andyV=17 



HOMOGCGC... /V=6and A/= 18 



LUMOGCGC... yV=5and/V=17 



LUMOGCGC... A/=6and/V=18 







FIG. 6: {\Afj,{t)\‘^) for HOMO and LUMO GCGC... if we initially place the carrier at the 1st monomer. [Left column] TV = 5 
and TV = 17. [Right column] TV = 6 and TV = 18. For TV even (|yl^(t)|^) are palindromes; for TV odd, this only holds for even 

M- While for holes |^| = |^| = 0.2, for electrons |^1 = |^| - 

*CG *CG 

magnitude for holes but rather similar for electrons. This leads to almost disrupted charge transfer if the number of repetition 
units is not integer i.e. for odd N. 


1.25, i.e., the hoping parameters^i^ are quite different in 
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LUMO ACAC... /V=5and/V=17 LUMO ACAC... /V=6and/V=18 
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FIG. 7: for HOMO and LUMO ACAC... if we initially place the carrier at the 1st monomer. [Left column] = 5 and 

N = 17. [Right column] N — Q and N = 18. We notice that for HOMO ACAC... where ]f^^] = ]f^’’ ], for N odd, (|A^(t)]^) are 
incidentally palindromes. This is a property mirrored from the stationary case since at this specific case ]u^fc]^ are palindromes. 
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Appendix A: matrix A 


A = 


■ j^bp(l) ,bp(l;2) 
^H/L ''ff/L 
,bp(2;l) pbp(2) 
^H/L ^H/L 


0 

bp(2;3) 

^H/L 


0 0 0 
0 0 0 


0 

0 


0 

0 


0 

0 


,bp(N-l-,N-2) pbp(Ar-l) ,bp(N-l-,N) 
^H/L ^H/L ^H/L 

„ bpiN-N-l) pbpiN) 

^ b-H/L ^H/L 


(Al) 


Appendix B: poly(dA)-poly(dT) as an example of type a' polymers 

In Fig. [10] we show some HOMO and LUMO properties of poly(dA)-poly(dT). We call Edge Group the first and the 
last monomer, and Middle Group the rest of the monomers. The total probability at the Edge Group, e{N) = 
and at the Middle Group, m{N) = [cf. Eqs. (1571) - (155)) ]. It seems that a power-law describes better the situation. 
This agrees with the claim that when every single hopping step occurs over the same distance, the hopping mechanism 
is described better by a power-law fi t^^i^^ . Since here we have the simplest periodic type one could imagine, the above 
claim evidently holds. In the last row we show u = kd versus d, which is less than 10^ m/s even for small oligomers. 
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FIG. 10: Hole (left column) and electron (right column) transfer in poly(dA)-poly(dT). [1st row] The total probability at the 
Edge Group, e{N) = and at the Middle Group, m{N) = [2nd row] For type a' polymers these probabilities are 

equally distributed among the monomers of these Groups. The probability at each of the members of the Groups is shown [cf. 
Eqs. (I37l) - (l38|l ]. [3rd row] The logarithm of the pure mean transfer rate k as. a. function of the distance from the first to the last 
monomer i.e. the charge transfer distance d = {N — l)x 3.4 A. [4th row] The logarithm of fc as a function of the logarithm of 
the number of monomers N. [5th row] The speed of charge transfer u = kd versus d. 
























